Spontaneous spin accumulation in singlet-triplet Josephson junctions 
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We study the Andreev bound states in a Josephson junction between a singlet and a triplet 
superconductors. Because of the mismatch in the spin symmetries of pairing, the energies of the 
spin up and down quasiparticles are generally different. This results in imbalance of spin populations 
and net spin accumulation at the junction in equilibrium. This effect can be detected using probes 
of local magnetic field, such as the scanning SQUID, Hall, and Kerr probes. It may help to identify 
potential triplet pairing in (TMTSF)2X, Sr2Ru04, and oxypnictides. 



PACS numbers: 74.50.+r effects 74.70.Pq 74.70.Kn 74.20. Rp 



Superconductivity with unconventional pairing, partic- 
ularly spin-triplet pairing, attracts a lot of interest in the 
condensed matter physics community and beyond. There 
is significant experimental evidence in favor of triplet 
pairing in the quasi-one-dimensional (QID) organic su- 
perconductors (TMTSF)2X P, i], ruthenate Sr2Ru04 
[31, and some heavy- fermion materials. In the recently 
discovered oxypnictide superconductors [J| , some experi- 
ments (observation of zero-bias conductance peak in tun- 
neling [5| and Hc2 exceeding the Pauli paramagnetic limit 
y) suggest a possible triplet pairing. Triplet [?[ and 
singlet |8| pairings were proposed in different theoretical 
models of oxypnictides. However, triplet pairing is not 
firmly established in any these materials. In this paper, 
we propose a new physical effect, which can provide use- 
ful information about spin symmetry of superconduct- 
ing pairing. We predict that electron spin accumulation 
should spontaneously develop at an interface between a 
singlet and a triplet superconductors. Other methods for 
detection of triplet pairing using superconducting junc- 
tions were proposed in Refs. [Sl, llO|. 

The predicted spin accumulation originates from the 
Andreev bound states at the interface between the singlet 
and triplet superconductors. Because of the mismatch 
between the spin symmetries of pairing, the spin up and 
down Andreev bound states have different energies. This 
results in different population of the spin up and down 
states and net spin accumulation at the interface. The 
preferred axis for the spin projection is determined by 
the vector d of the triplet pairing, and the sign of the ac- 
cumulated spin is determined by spontaneous symmetry 
breaking. The resulting magnetization can be detected 
using local probes of magnetic field, such as the scanning 
SQUID, Hall, or Kerr probes. If a voltage V is applied 
to the junction, then the magnetization would oscillate 
with the Josephson frequency 2eV/h. In this paper, we 
present calculations in two cases: for the non-chiral p^- 
wave pairing, relevant to the (TMTSF)2X materials [llj . 
and for the chiral px + ipy pairing, relevant to Sr2Ru04 
[12] ■ For oxypnictides, the calculations would be more 



complicated because of the multiple bands [8|, but the 
result should be qualitatively the same. Although the 
Andreev bound states [13| in the singlet-triplet Joseph- 
son junctions [IJI were studied in literature before, the 
spin accumulation effect was not recognized, except in 
Ref. 15] for the special case of equal energy gaps. 

Let us consider a Josephson junction between an s- 
wave and a p-wave superconductors located at a; < and 
X > respectively, as shown in Fig. [T] The two supercon- 
ductors are separated by a narrow insulating barrier at 
X = Q, which is modeled by the delta-function potential 
U{x) = Uo6{x). The interface between the superconduc- 
tors is assumed to be smooth, so that the electron mo- 
mentum parallel to the interface fcii = {ky,kz) is a good 
quantum number. The interface plane is perpendicular 
to the planes of Sr2Ru04 or chains of (TMTSF)2X, as 
shown in Fig. [1] The pairing potential of the singlet s- 
wave superconductor on the left (L) side is 



{ck,aC-k,s) fx Aisgn(o-) = A^ 



(1) 



where Ck,cr is the destruction operator of an electron with 
the momentum k and spin a =t, J,. Here, sgn(cr) = -l-(— ) 
and a =| (j) for a =t (j). The pairing potential of the 
triplet p-wave superconductor on the right (R) side is 

{Ck.aC-k,a') CC A2 i(7,,(o- ■ d) f (k) e'^" . (2) 

Here <t are the Pauli matrices in the spin space, d is a 
unit vector characterizing the spin polarization of the 
triplet superconductor, and (j>o is the U(l) phase dif- 
ference across the junction. The function /(fc) rep- 
resents the orbital symmetry of the pairing potential: 
/(fe) = {kx + iky)/kF for the chiral Px + ipy pairing and 
/(fe) = kx/kp for the non-chiral Px pairing, where kp is 
the Fermi momentum. We assume that the vector d has 
a uniform orientation independent of fe. By selecting the 
spin quantization axis z along d, Eq. ([2|) is simplified as 



(cfc^.c_fc,^)«A2/(fc)e'^«EEA^. 



(3) 



In this representation, both singlet ^ and triplet ([3|) 
pairing potentials couple electrons with opposite spins. 
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FIG. 1: A Josephson junction between a single and a triplet 
superconductors. The thin black solid lines represent chains 
for (TMTSF)2X or planes for Sr2Ru04. The blue dotted 
lines indicates the localization lengths of the Andreev bound 
states. The cone represents a probe of the local magnetic field 
produced by spin accumulation at the interface. 



Electron states in a superconductor are described by 
the Bogoliubov operators 7, which are related to the elec- 
tron operators c by the following equations [l6| 



7<Tnfc|i = / fia;[<„fe||(x)c^fc||(a;) +<„fe||(x)c^^^^(a;)], (4) 
Cak^fix) - ^K„fc||(a;)7^„fc|| + ^fffe,, (a;) 7^_^ ], (5) 



where fe|| — — fc||, and n is the quantum number of the 
Bogoliubov eigenstates. The two-component wave func- 
tions Vo-n(a:,fc||) = [uankiiix),Vankiiix)] are the eigen- 
states of the Bogoliubov-de Gennes (BdG) equation with 
the eigenenergies Ea-nk« 



Ho + U{x) A^(a;,fc^,fc||) 
A;(x,fc„/C||) -Ho-Uix) 



tpaix) = E„%l^a{x), 



(6) 
where we omitted the labels n and feii for -0 and E to 

shorten notation. Here Ho = £{kx, fc||) ~ M is the electron 
dispersion relation with kx — ~idx and the chemical po- 
tential (U, and we set fi = \. In Sr2Ru04, the main Fermi 
surface is circular, so we take e = {k^ + ky)/2m, where 
m is the effective mass. In QID conductors, the Fermi 
surface consists of two open sheets perpendicular to the 
chains, so we take e — fc^/2m — if, cos{bky) — tc cos(cfcz), 
where the first term represents motion along the chains, 
and tt, tc, b, and c are the interchain tunneling ampli- 
tudes and spacings [17|. Notice that the spin projection a 
is a good quantum number for the Bogoliubov quasipar- 
ticles Q, and the BdG equations © separate for a =t 



and |. The pairing potential Ao- in Eq. ([6]) is given by 
Eq. (P) for X < and by Eq. ^ ior x > 0. 

The wave functions ip^ and ipa '^^ the left and right 
sides of the junction satisfy the standard boundary con- 
dition at a: = obtained by integrating Eq. ^ over x 
from -0 to +0 

^^(0) = ^^(0), {dx^/j^ - dx^!;)x=o = 2mUo^!;{0). (7) 

We are interested in the subgap bound states of Eq. ([6]) 
with energies \E\ < Min[Ai, A2], which are localized near 
the junction. Such localized solutions can be obtained as 
a superposition of the wavefunctions for the right and left 
moving quasiparticles [la |: 



0(^\ - 



rAx) 



A^ 



^ikpx 






-B 



/3 -ikpx 



.fi 



(8) 

Here the superscript fi = R,L labels the wave functions 
on the right and left sides of the junction. The sub- 
script a = ± in Uaa or v^a denotes the right (a = -I-) 
or left (a = — ) moving quasiparticles. The parameters 



R{L) 



V 



A2 

^2(1) 



E'^/vp determine the inverse lo- 
calization lengths of the bound states inside the right 
and left superconductors. The variables kp and vp are 
the X components of the Fermi momentum and Fermi 
velocity, which, generally, depend on fc||. For Sr2Ru04, 

= kp/m. 



— \i kp 



fc2 and Vp 



For QID conductors, 



kp = kp+2tbCos{bky)/vp+2tcCOs{ckz)/vp and Ui? « vp, 
where vp — kp/m. The coefficients m^„ and v^^ are de- 
termined by substituting the right and left moving terms 
into Eq. ^ away from the junction. They satisfy 



'laa 



m'3 



E„ ~isgn{a)vpK^ 

Ag 



(9) 



where A^ are given by Eqs. (jT]) and ^. Notice that all 
variables in Eqs. ([8]) and ([9]), generally, depend on fc||. 

Substituting Eq. ([8]) into Eq. ([7]) , we obtain a set of 4 
linear homogeneous equations for the coefficients A^ and 
B^ . A condition for non-zero solutions requires vanishing 
of the determinant of the corresponding 4x4 matrix, 
which yields the following equation 






vi-] 



1 -^(fcll 



(10) 



Here D(/c||) = 4/(4 + Z'^{k\\) is the transmission coef- 
ficient of the barrier. For QID conductors, Z{k\\) ~ 
Zo = 2Uo/vp is independent of the transverse momen- 
tum. For Sr2Ru04, Z{k\i) = Zakp/Jkj, - fc^. Substi- 
tuting Eq. ([9]) into Eq. pO|) . we obtain an equation for 
the energies of the Andreev bound states 



^sin($fc ) -6cos($fc ) =sgn(CT) sin{(l3o) D{k»). (11) 




FIG. 2: (a) The energies Ef (dashed lines) and Ei (sohd lines) 
of the Andreev bound states vs. the phase (po between the s 
and p^-wave superconductors {D = 0.8 and A2 = 0.5Ai). (b) 
Magnetic moment m^ per unit area of the interface, Eq. (|12|l . 
vs. (j>o for T/A2 — 0.01 (solid line) and 1 (dashed line). 
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FIG. 3: The energies E-\{ky) (blue dotted lines) and Ei{ky) 
(red dash-dotted lines) vs. ky for the phase (po — 7r/2 between 
the s and p^ + ipy superconductors {Zo = 1). For (/>o = and 
TT, E^(ky) = Eyiky) is shown by the black solid line. 



Here (f){k) is the phase of the function /(fe) in Eq. ([3]), and 
*&fe|] — [4>{kF, k\\) — 0(— ^F7fe||)]/2 is a half of the phase 
difference between the points on the Fermi surface con- 
nected by specular reflection from the barrier, selected so 
that < <I> < TT [l2|. The coefficients A and B are 



A 



[2 ~ i?(fc||)] e.i a/1 - €^2 + I?(fc||) e,2 J I - e^^. 



B = -D{k\\)e,,e,^ + [2-D{k^)]J{l-el,){\-el^ 



where £0-1(2) — Ea-{kf, 0o)/^i(2) < f is the dimensionless 
energy of the bound state. In the case Ai — A2, the 
equations simplify and reproduce the results of Ref. [l^ . 

Notice that E^ depends on the spin index a only 
through the right-hand side of Eq. (fTTj) , which is invari- 
ant under the transformation a ^> a and 0o -^ 00+"', so 
Ea{k\\\(t)) = Ea{k\\]4) + tt). This relation can be under- 
stood by noting that the singlet pairing potential ^ has 
opposite signs for u =] and |, whereas the triplet pairing 
potential ^ has the same sign. So, a phase difference 00 
across the junction for spin-up quasiparticles implies the 
effective phase difference 0o + tt for spin-down quasiparti- 
cles, which explains the above-mentioned invariance. We 
see that the bound states energies E^ and Ei are gener- 
ally different for a given (pQ, which is the key point for 
understanding of spin accumulation [f8!|. 

For thepa;-wave pairing, we have f{—kp) = —f{kp), so 
$fc|i = 7''/2, and Eq. ([TT|) reduces to ^ = Dsgn(cr) sin(0o), 
where D = A/{A + Zq). In this case, the energies E^ 
are independent of fc|| and are plotted vs. 0o hi Fig. [2^. 
Depending on 0o, there are two or four of such states 
for each A;||. The spin up and down bound states have 
opposite energies i?| = — £'|, as shown by the dashed and 
solid lines in Fig. [2^. The difference between the Fermi 
populations of the spin up and down states gives a net 



magnetic moment niz per unit area of the interface [20| 

(12) 



^^tanh^'^^^'^") 



2hc 



2T 



where /is is the Bohr magneton, T is the temperature, 
and the prefactor 1/2 compensates for double-counting 
lllUla|. A plot of TTiz vs. 00 is shown in Fig.[2l3 for two 
different temperatures. At low T, the magnetic moment 
is close to [IbI2 per chain, since only the lower energy 
state is populated. In an open circuit, the value of the 
phase 00 is determined by minimization of the total en- 
ergy of the system. For the energy levels shown in Fig. [2^, 
the minimum is achieved at either 0o = 7r/2 or 0o = 37r/2 
(the same as 0o = — 7''/2) [19]. The system spontaneously 
breaks the symmetry and selects one of the two energy 
minima with negative or positive magnetization. 

For the chiral j)x + ipy pairing, we have 0(fc) = 
arctan(fcj,/fc2;) and <I>fc = 0(fca) + '^ I'^i so that sin(f>fe = 
fcxI/^F and cos$fc = —kyjkp. Then, Eq. (fTTj) gives the 
energies Ea-{ky; 0o) dependent on the transverse momen- 
tum ky because of the broken time-reversal symmetry. 
The plot of Ecy{ky] 0o) vs. ky for several values of 0o is 
shown in Fig. [31 The energy splitting between E-^ and 
E^ is maximal at 0o — ±7r/2 and vanishes at 0o = and 
TT. The imbalance between the spin up and down popu- 
lations produces the net magnetic moment rriz per unit 
area of the interface 



2c 



Y^ sgn(cr) 



"2^ 



np 



Ea{ky](j)Q) 

T 



(13) 



where np is the Fermi distribution function, and c is the 
interplane distance for Sr2Ru04 [20|. The plot of niz 
vs. 00 is shown in Fig. [5] for two different temperatures. 
The minimum of energy is achieved at 0o = ±7r/2, so the 




FIG. 4: Normalized magnetic moment tti^ per unit area of 
the interface, Eq. (|13|) . vs. <jf>o for T — O.OIA2 (solid line) and 
T = A2 (dashed line) for a junction between the s andpx+ipy 
superconductors (p = kF/2nc). 



system spontaneous breaks the symmetry and selects one 
of the two optimal values for (Jjq . 

Figs. [2Jd and S] show that the magnetic moment tti^ 
changes sign when (J>q crosses tt. If a bias voltage V is 
applied across the junction, it would make the phase dif- 
ference time-dependent: (IJoit) = 'ieVt. Then, the mag- 
netization ruz at the interface would oscillate with the 
time period ir/eV . We assume that the oscillations are 
slow enough for the spin population to remain close to 
the thermal equilibrium at each moment of time. 

Finally, we discuss possible experiments for detection 
of the spontaneous spin polarization. A schematic ex- 
perimental setup is shown in Fig. [1] Assuming semi- 
infinite geometry in the z direction, the magnetization 
at the junction and the magnetic field can be estimated 
as i? ~ /io^_B/OK/4, where k ~ A2/VF is the inverse pen- 
etration depth of the bound states in the case Ai :§> A2, 
p = I /be for (TMTSF)2X, and p = kpI'l'KC for Sr2Ru04. 
This magnetic field can be measured using a scanning 
SQUID or Hall microscope [U [ll]. For (TMTSF)2X, 



we have h = 0.77 nm, c = 1.35 nm and 



0.6 /im, 



which gives B ~ 0.3 G and a magnetic flux 0.06$o 
(where $0 = hc/2e) through a square scanning SQUID 
loop of the size I = 10 /im. For Sr2Ru04, we have 
kp = 7.5 X 10^^ m^^, K~^ = 66 nm and c — 1.3 nm, 
which gives the field 0.7 G and the flux 0.15$o- The es- 
timated magnetic fields are well above the typical Hall- 
probe sensitivity of 80 mG at 1 Hz 2l|. However, chiral 
superconductors are also expected to have an additional 
magnetic field due to the charge currents carried by the 
chiral Andreev bound states [23|. A Josephson junction 
between Auo.slng.s and Sr2Ru04 was scanned using the 
SQUID and Hall probes in Ref . [221 , but no spontaneous 
magnetic field was detected. A reason for the negative ex- 
perimental result remains an open question. Local mag- 
netization can be also detected optically using the Kerr 



angle rotation 24|. This effect was observed in Ref. 25 1 
in the bulk of Sr2Ru04 due to the orbital time-reversal 
symmetry breaking in the Px + ipy state. However, the 
Kerr experiment has not been performed in a scanning 
mode at a junction with a singlet superconductor. 

In conclusion, we have shown that a Josephson junc- 
tion between a singlet and a triplet superconductors 
should exhibit spontaneous spin accumulation due to 
mismatch of the spin pairing symmetries. The vector 
of the accumulated spin points along the vector d of the 
triplet superconductor. The sign and magnitude of the 
spin depend on the phase difference (f>o between the su- 
perconductors. In equilibrium, the system spontaneously 
breaks symmetry and selects one of the two values of (po 
that minimize total energy and maximize spin accumula- 
tion. When a bias voltage V is apphed to the junctions, 
the accumulated spin oscillates in time. The magnetic 
field produced by the accumulated spin can be detected 
using the SQUID, Hall, or Kerr local probes. 
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